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Abstract 

We investigate the propagation of electromagnetic waves in finite photonic 
band gap structures. We analyze the phenomenon of conduction and forbid- 
den bands and we show that two regimes are to be distinguished with respect 
to the existence of a strong field near the interfaces. We precise the domain 
for which an effective medium theory is sounded. 

I. INTRODUCTION 

The theoretical and numerical study of photonic band gap materials (see |I[ for an 
exhaustive bibliography) may be led both from the point of view of the spectrum, with the 
computation of the band structure by using Bloch wave theory and the related quantities 
(Density of States and so on), or from the point of view of the scattering theory ||. In 
case of a scattering experiment, which of course involves a finite structure, Bloch waves 
are not a priori sufficient to describe the electromagnetic field because of the existence of 
evanescent and anti-evanescent waves, linked to the boundary of the device. Of course 
these non-propagating modes are of prime importance in the band gaps, as they represent 
entirely the total field inside the photonic crystal, but they may also induce strong effects 
for frequencies lying inside the conduction bands. 

In the context of electromagnetic optics and the spectacular effects that can be obtained 
using photonic crystals (the so-called ultra-refraction effects [[^[|],|5],|C|[7||| ) a theoretical ap- 
proach only involving equivalent medium theories, group velocity |||4]|| and more generally 
quantities only derived from the band structure is certainly incomplete. Previous studies 
have shown the importance of considering the isoenergy diagram and not only the disper- 
sion diagram . In the present work, our aim is to analyse the relative importance of the 
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evanescent waves (the near-field) for describing the electromagnetic field inside a photonic 
crystal. More precisely, we study the ratios of the projection of the field on the propagating 
and non-propagating modes, and we precise the conditions under which Bloch waves can 
describe entirely the scattering behaviour of a photonic crystal. 



II. PROPAGATING AND NON-PROPAGATING MODES 

The photonic crystal under study is made of a stack of gratings. For the numerical 
applications we use square dielectric rods in a dielectric matrix (fig. 1). The basic layer of 
the device is made of one grating, consisting of a slab of rods between two homogeneous 
slabs above and below (fig. 1): it covers a band Ix [—h/2, h/2]. The period of the grating 
is d. We use time-harmonic fields with a time dependence of e~ luJt and k — ^ denotes 
the wavenumber in vacuum. We have chosen this particular geometry (square rods and 
rectangular symmetry) and polarization for the sake of simplicity (both on the numerical 
level and for the theoretical exposition). However the theoretical framework is quite general 
and works as well for more complicated situations (see |TIJ for an exposition of the transfer 
matrix method for 3D obstacles). 

Our aim is to characterize the field inside the photonic crystal. Due to the transla- 
tional invariance of the problem, we look for pseudo-periodic fields in the x direction: 
u(x + d,y) = exp (iad) u (x, y), where a belongs to the interval Y = [—ir/d,7r/d[. In a 
diffraction problem, where the grating is illuminated by a plane wave under the incidence 8 
(fig. 1), we have a = ksin.8 mod(7r/cf). From grating theory ||11|| , we can define the trans- 
fer matrix of the basic layer, which is the operator T( a ,\) linking (u(x, h/2), d y u (x, h/2)) to 
(u (x, —h/2) , d y u (x, —h/2)) (in ordinary differential equations theory this is the monodromy 
matrix). There are many very good numerical methods for computing the field diffracted 
by a grating and hence the transfer matrix, even for stacks of gratings |12| , |i3| , |14|1 . 

On the upper (resp. lower) side of a basic cell, we expand the field on a Rayleigh basis: 

u(x, h/2) = J2( A n+ K) e tanX , u(x, -h/2) = £ (S+ + B~) e fa "*. (1) 

n n 

where a n = a + nK, K = 2£. The values of the normal derivatives d y u\ y=±h ^ 2 write: 
d y u(x, h/2) = £ i/3 n (A+ - A-) e ia - x , d y u(x, -h/2) = £ ifl n (B+ - B~) e ia " x . 



The knowledge of A ± = {A^} n (resp. B ± = {B^}) gives the value of the derivatives. 
Therefore, rather than computing the monodromy matrix as defined above, we compute the 
matrix T( a ,\) such that 

(2) 

The point is to analyze the spectrum of T. For symmetry reasons, the spectrum 
sp (T) of T is invariant under r — ► r _1 (this is easily seen in case of a lamellar grat- 
ing, but the proof is slightly more involved in case of a y-dependent medium), then we 
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can distinguish between eigenvalues of modulus 1 that are necessarily finitely many, and 
eigenvalues that do not belong to the unit circle of the complex plane. Let us denote 
e ll3h ,P G — ?, \ , an eigenvalue of T( a ,\) of modulus one, and ip an associated eigenvector 

then we have T ,ip~^J = e t/3h V' - ) or e l se: + = e l ^ h ip (x,y) this means 

that (x,y) — * e %ax ip (x,y) is a Bloch wave associated to the Bloch vector (a, (3). That way 
we can easily compute the dispersion curves at a given wavenumber k. Moreover, we can 
also compute the non-propagating modes inside the crystal: they correspond to eigenvectors 
associated with eigenvalues that are not of modulus one. We have thus obtained a decompo- 
sition of the modes by means of a family of monodromy operators parametrized by a G Y . 
As it as been said before, from the scattering point of view, the parameter a is equal to 
k sin 9 so that we study the spectrum of T(e,\) — T( k sme,\)- 
We can now give the following definitions. 

Definition 1 We call relative gap an interval of wavelengths Ig, at a given incidence 9, 
for which T^a) has no eigenvalues of modulus one and we call relative conduction band 

an interval Bg of wavelengths where T(# jA ) does have eigenvalues of modulus 1. 



Definition 2 A total gap corresponds to the intersection of incident dependent gaps ( and 
may be void). 

III. ANALYSIS OF THE SPECTRUM 

As it has already been stated, at a given wavelength, Bloch waves are not sufficient to 
compute the scattering properties of the crystal. In order to quantify the relative importance 
of the evanescent waves, we need to be able to get a decomposition of the field. 

A. Decomposition of the field 

Once the electromagnetic field is known on the upper face of the crystal (through the 
coefficients A ± , obtained from a rigorous numerical method), it is possible to expand it 
on the various modes that exist in the grating layer. More precisely, except on a set of 
wavelengths of zero Lebesgue measure, matrix T( a ^) can be put in diagonal form: 

T( a ,A) = T p ©T e ©T a 

where T p is a finite rank operator corresponding to propagative waves and T e ,T a correspond 
to the evanescent and anti-evanescent modes. 

With this decomposition, the vector ip = ^ + ,-^~^ writes ip = ipp^ipe^^a- Whence we 
define the branching ratios n p (resp. 7r e , n a ) of the field on the propagating (resp. evanescent, 
anti-evanescent) modes by: 

p ~ NW e ~ a ~ N(il>) [ } 
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where iV(^) = ||^|| 2 + ||^e|| 2 + ||^ ( 




The point of the above decomposition is to quantify the relative importance of the various 
modes in the total field existing in the crystal, in order to understand to what extend the 
field is not solely described by Bloch waves. 

B. Cut wavelengths and classification of the conduction bands 

Let us now turn to some numerical computations. The relative permittivity of the rods 
is 62 = 9 and e% = e ext = 1, the geometric parameters are h = 2.8, d = 2.8, hi = 1.9, d\ = 1. 
The structure is made of one basic layer and we choose 9 = 30° and s-polarized waves. In 
fig.2 (a), we give the absolute values of the eigenvalues of T e (A) versus the wavelength. 
The conduction bands are the regions with a horizontal straight line (\fi\ = 1). For each 
wavelength A there is a finite, possibly empty, set of eigenvalues of modulus one |e^ n( - A ^| 
and an infinite set of eigenvalues that do not belong to U = {z G C, \z\ = 1} (in fig. 2 (b) we 
have plotted the real part of the spectrum of T^^, where the finite number of propagating 
modes may be observed). We have also plotted in fig.2 (c) the dispersion diagram (/3, X/d). 
The comparison with fig. 2 (a) shows that the consideration of the complete spectrum of 
T(e,\) , i-e. with the non-propagating modes, allows to understand that one should not 
treat on a different foot propagating and non-propagating modes because they are really the 
same physical entities, behaving differently according to the wavelength. When 9 varies in 
]— 7r/2, 7r/2[, the local gaps vary as shown in fig. 3. 

This first example is rather generic and shows that within a given conduction band, hence 
locally, it is possible to define continuous sections A — > /i n (A) G sp (T^^ representing the 
evolution of the eigenvalues of the monodromy operator with respect to the wavelength. At 
some values of the wavelength however, these sections may encounter a bifurcation, or cut- 
off: the eigenvalue leaves U and the associated modes give rise to an evanescent mode and an 
anti-evanescent mode. At such a branch point, the section A — > fi n (A) is not differentiable 
and may cross other sections. As a consequence, a global description of the sections is not 
possible in that case. This problem is quite a complicated one, for which there is a general 
theory pO| , pT| , but even deriving a specific theory for our particular situation is quite a big 
task and beyond the scope of this work. 

However, an easy simplification can be obtained by noting that the set of eigenvalues 
being invariant under r — > r _1 , it is natural to consider the quotient space sp (t^^A / ~ 
for the equivalence relation fii ~ fi 2 if ^1/^2 = 1 ; which amounts to identify two eigenvalues 
that are inverse one of the other. This operation gives a nice simplification but still does 
not allow to define global sections (a more detailed account of this situation will be given 
elsewhere ||15||). 

Remark 3 A very simple example of branch point is the extinction of a diffracted order in 
grating theory. Another elementary situation is that of a stratified medium ( a Bragg mirror 
for instance) in which case there are only two propagative modes in the conduction bands and 
one evanescent and one anti- evanescent mode inside the gaps. A realization of the quotient 
space is obtained by considering ||tr (T^x)j , A G 1R + |. In that case, this very set defines 

a global section and the quotient space sp (T^ x)) / ~ is a trivial fibred bundle. 



4 



For a given incidence, a gap is then an interval of wavelengths over which all the prop- 
agative eigenvalues have encountered a bifurcation. In fig. 2 (a), we have this situation in 
the interval (1.32, 1.42). 

When the wavelength tends to infinity, it is known that the device finally behaves as 



a homogeneous slab [|T^,|T^,|T^,|T9|] , and then there are only two propagative modes (up and 
down), which means that all "sections" finally bifurcate definitely (see fig. 2 (a-b) for 
X/d > 2.22), except the one corresponding to the homogenization regime. In that case there 
are still evanescent (and anti-evanescent waves) but with a very huge damping exponent so 
that 7r e and ir a are small. 

However, before that regime, eigenvalues may experience a local bifurcation: that is they 
leave U over a finite interval but finally come back on it (in fig. 2 (b) this situation happens 
over the interval (1.63, 1.87)). What is important to note is that such a local bifurcation may 
affect only one eigenvalue so that, whithin a conduction band there may be evanescent field 
coming from such a bifurcation, hence with a small damping exponent (this happens over 
the interval (1.32, 1.42) in fig. 2 (a-b)). This leads us to distinguish between both regimes 
and give the following definitions. 

Definition 4 A conduction band is said local if among the evanescent modes within this 
band there is at least one mode corresponding to a local bifurcation. A conduction band is 
global if it is not local. 

Let us now give some numerical examples of the various regimes described above. We 
give in fig.4 the absolute values of the eigenvalues of matrix TVa) (normal incidence) for a 
one layer structure with e ext = 2.26, e± = 1, Si = 4 h = 1, d = 1, di — 0.5 and the projection 
ratii 7r Pje ,a- The region (1.28, 1.37) corresponds to a local conduction band, i.e. in which 
there is a local bifurcation of an eigenvalue. We see that the part of the field on the non- 
propagating modes is not at all negligible so that the field cannot be described solely by 
Bloch modes. On the contrary, for the interval (1.72, 1.8) the conduction band corresponds 
to a global bifurcation of two eigenvalues of modulus one and in that case the damping 
exponents are exponentially growing, so that almost Bloch waves only contribute to the 
description of the field. 

A natural question is to know to what extend these situations persist when the number of 
layers is increased. We have computed the values of the branching ratios when the number 
N of layers is N = 2, 4, 6. The results are given in figures 5 (a-b-c). It can be seen that the 
branching ratios ir e and n p exhibit an oscillatory behavior with respect to A, which is probably 
linked to the excitation of resonances, but that the fraction of electromagnetic energy that 
is carried by the evanescent waves is not diminished. This means that a non negligible part 
of the field is localized near the interfaces, which can have substantial consequences on the 
propagation of a beam inside the structure (this situation will be analyzed in a forthcoming 



paper |17|) but also on the local density of states for photons. 

As an application, suppose now that the wavelength belongs to a global conduction 
band and is such that there is only one diffracted and one transmitted order; then from the 
knowledge of the transmission and reflection coefficients, it is possible to compute the (2 x 2) 
monodromy matrix of the device. In that case, we can obtain the superior envelope of the 
transmitted energy by considering only the transfer matrix of one basic layer [E3] . We have 
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shown in another article p3] that in a layer characterized by a 2 x 2 transfer matrix, the 



reflected and transmitted coefficients for N layers can be obtained in close form : 

where /i is an eigenvalue of T p = (t^) associated with an eigenvector u = (■u 1 ,-u 2 ), an 
eigenvector associated to is denoted by v = (t>i,t>2) and, denoting q(xi,x 2 ) = Ip^l+xl > 

/5o = — cos(0), functions / and g are defined by 
A 

if (M)6G g(.k,0) = q(y), f(k,6)=q(w) 

if r*. m *= r I * ^ = g (v) ' ; ^ = g (w) if |g (v)l < |g (w)l (5) 

1 ' J U(M) = g(w),/(M) = ff(v) if [g(w)| < |g(v)| 
where G = {(&, 6 1 ) , |tr(T p )| > 2} and B = {(k, 6) , |tr(T p )| < 2} . 

The superior envelope i?^ of the reflected energy, and conversely the inferior envelope of 
the transmitted energy are given by |22] : 



_ 4-tr(T p ) 2 

J- nc 



(ti2/9o - tsi^" 1 ) (6) 



- oo 



A direct application of these formulas show a very accurate result in fig. 6 (a-b) for a global 
conduction band. 

We can conclude by the following: 

Proposition 5 Within a global conduction band, the field inside the crystal can be repre- 
sented by Bloch waves only. 



IV. CONCLUSION 

We have shown that it is important to distinguish between various kinds of conduction 
bands: there may be non propagative modes that result from the local bifurcation of a 
propagative mode or all the non propagative modes may be made out of global bifurcation 
of propagative modes. In the first case, an important part of the field inside the structure is 
made of non-propagative modes, but in the second case, the field writes in terms of Bloch 
waves only. Especially near a band edge, one should be very careful before deriving the 
behavior of the field solely by looking at the dispersion diagram: it does not take into 
account the evanescent waves. These results might be useful in studying beam propagation 
and superprism effects in photonic crystals [p4j] and also the phenomenon of spontaneous 



emission in a finite dielectric structure, for in that case the density of modes for photons 
derived without precaution from Bloch theory is certainly false because the atom modes can 
couple to non-propagating radiation modes |p5| . 
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Figures captions 

Figure 1: Sketch of the photonic crystal. 

Figure 2: (a) Absolute value of the eigenvalues of Tg (A) for e\ = 9, Si = e ext = 1, h/d — 
1, hi/d = 0.68, d\/d = 0.35, s-polarized waves. 

(b) Real part of the spectrum of Tg (A). The arrows indicate a local and a 
global bifurcation. 

(b) Dispersion diagram (/?, X/d) for the parameters of figure 2 (a). 
Figure 3: Evolution of the gaps with 9. The white regions correspond to band gaps. 
Figure 4: Branching ratios for e ext = 2.26, 62 — 4, e\ = 1, h/d — 1, d\jd = 0.5, s-polarized 
waves. 

□ : propagating ratio (ir p ) , A: evanescent ratio (7r e ) , *: anti-evanescent ratio 

Figure 5: Spectrum and branching ratios for the parameters of figure 4. 

(a) iV = 2 (b) N = 4 (c) N = 6. 
Figure 6: Reflected energy and its enveloppe given by @, with iV = 20. 

(a) E\ = 9, Ei = Eext = 1, h/d = 1, hi/d = 0.68, d\/d = 0.35, s-polarized waves. 

(b) E ext = 2.26, £2 — 4, £1 = 1, h/d = 1, d\/d = 0.5, s-polarized waves. 
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